Optical Dichroism by Nonlinear Excitations in Graphene Nanoribbons 
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The honeycomb carbon structure of graphene and nanotubes has a dynamics which can give 
rise to a spectrum. This can be excited via the interaction with an external electromagnetic field. 
In this work, non-linear waves on graphene and nanotubes associated with the carbon structure 
are investigated using a gauge model. Typical energies are estimated and there scaling with the 
nanoribbon width investigated. Furthermore, the soliton-photon interaction depends on the incident 
photon polarization. In particular, we find that the nanoribbon is transparent when the polarization 
is along the largest length. Relying on the scaling with the width, we suggest a way to experimentally 
identify the soliton waves in nanoribbons. 

PACS numbers: 71.10.-w,72.80.Vp,11.10.Kk 
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I. INTRODUCTION AND MOTIVATION 

The physics of nano-materials has been the subject 
of intense investigation in the past years. In particular, 
graphene and nanotubes^"— have been studied due to its 
potential technological applications and its peculiar elec- 
tronic properties, see^ii. 

A popular approach to understand graphene and nan- 
otubes is to use tight-binding or Hubbard like models to 
describe the behavior of electrons and holes. Despite its 
success, these type of models do not take into account 
the dynamics of the honeycomb array of carbon atoms. 
Therefore, dynamical effects connected directly with the 
carbon structure are not taken into account by those type 
of models. 

A chiral gauge model for graphene was suggested in^. 
In^ it was studied its generalisation to include other 
gauge groups, not necessarily chiral invariant. The mod- 
els are able to accommodate the main features of the 
electronic properties in graphene, like the creation of 
mass gaps or quantum Hall effects, and parameterizes 
the dynamics of the carbon background structure via a 
charged scalar field (p and spin one fields A^^ . The model 
allows to go beyond the electronic properties and to in- 
vestigate phenomena which arc related with the dynam- 
ics of the carbon crystal structure. Indeed, the authors 
explored the non-linear dynamics of (p and to com- 
pute mass gaps, to derive finite energy vortex solutions, 
which can lead to flux quantisation and Aharonov-Bohm 
effects (see e.gJ^ and references therein) and to check 
that the gauge theory is compatible with the experimen- 
tally observed electronic quantum Hall effect. Further- 
more, the fermionic zero modes associated to the finite 
energy vortex configurations^ are solutions of a one di- 
mensional Schrodinger-like equation with a scale invari- 
ant 1/z^ potentiaiii. Certainly, the theory can accommo- 
date other phenomena not explored so far and, in par- 
ticular, phenomena whose dynamics is linked with the 
honeycomb array of carbon atoms. 



The complex scalar field (p and model the dynam- 
ics of the carbon crystal structure, on top of which elec- 
trons and holes live. Certainly, the electronic properties 
of graphene are sensitive to the (p and A^ configurations. 
Changing the configuration of the bosonic fields, the elec- 
trical and thermal transports properties associated with 
the fermionic degrees of freedom change accordingly. For 
example, phonons are associated with fluctuations of p 
and Af^ around a given conflguration and contribute to 
the graphene specific heat. Then, the investigation of the 
thermal properties of graphene can help understanding 
the bosonic part of the gauge model. On the other hand, 
the scalar field p being a charged field, it couples directly 
to the electromagnetic field and can contribute to the 
graphene optical properties. Therefore, within the gauge 
model o:^— , it is of paramount importance to understand 
how the bosonic degrees of freedom can be excited via the 
coupling with an external field. If one is able to control 
the p and A^^ configurations by coupling with external 
fields, one can achieve a better control of graphene and, 
in general, nano-materials properties. 

From the pure field theory point of view, p couples 
directly to the fermions and, for example, can trigger the 
generation of a fermionic mass gap. Given that p is a 
dynamical field, the model allows us to define regions 
where the mass gap vanishes, while in other regions the 
fermions acquire a finite mass with the corresponding 
decreasing of its mobility, i.e. from the point of view 
of the gauge model the mass gap is, itself, a dynamical 
quantity. One would like to understand what type of p 
configurations are allowed and how they impact on the 
electronic properties. 

The self-couplings of p allow for solitonic solutions of 
the fields equations. In this work, we will use the term 
solitonic solution to describe any solution of the non- 
linear differential equation of motion for p. The field fluc- 
tuations around the soliton can give rise to phonons with 
new types of dispersion relations. One expects contribu- 
tions to the optical and thermal properties of graphene 
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coming from tlic solitonic waves. Furtlicrmore, lias a 
non-vanisliing electric charge and the non-linear waves 
couples with an external electromagnetic field. This cou- 
pling opens the possibility of generating these non-linear 
modes by exciting graphene with a laser or a properly 
tuned electromagnetic field. 

The gauge model suggested in^ is an effective quan- 
tum field theory for the dynamics of a 2-|-l dimensional 
system. In principle, it can be applied to any 2D nano- 
materials, other than graphene. In this work we will 
also consider nanotubes. From the point of view of the 
gauge theory, the main difference between the two sys- 
tems is the geometry of the two carbon based materials. 
If graphene is associated with a 2D flat sheet, a nanotube 
is a two dimensional material having the same type of 
unit cell but with a cylindrical geometry. 

In this paper, we investigate the classical solutions of 
the non-linear bosonic equations of motion of the gauge 
model? for graphene and nanotubes. Besides computing 
explicitly no n- linear configurations for ip and A^, we also 
discuss how these nonlinear modes couple with an exter- 
nal electromagnetic field. In particular, we find that the 
coupling between the nonlinear modes and the photon 
field displays dichroism, i.e. the ability to absorb light 
depends on the polarization state of the incoming photon. 
Dichroism with a strong absorption anisotropy from the 
electronic intraband transitions requiring a finite chemi- 
cal potential was investigated inii. For these two cases, 
the pattern of the coupling with the polarization of the 
photon is similar, i.e. the graphene nanoribbon is trans- 
parent when the polarization is along the largest length 
of the carbon material, but also the energy range where 
dichroism is expected to be seen is in the far infrared 
and terahertz frequencies. If the intraband transition re- 
quires a non-vanishing chemical potential to generate 
dichroism, the coupling with the nonlinear waves is inde- 
pendent of [I. This observation enables to experimentally 
distinguish between the two mechanisms. 

The paper is organized as follows. In section |ll] wc 
review the gauge model. The coupling with an exter- 
nal electromagnetic field is discussed, the equations of 
motion and the Hamiltonian for the gauge model are de- 
rived. In section Hm the solitonic solutions of the classical 
equations are discussed in the present context. In section 
II VI we investigate the soliton-photon transition amplitude 
and decay rate. Finally, in section |V] we resume and con- 
clude. 



II. GAUGE MODEL FOR GRAPHENE AND 
NANOTUBES 

For completeness, in this section we resume the main 
features of the gauge model described in^. The inter- 
ested reader should look at the cited paper for a detailed 
discussion. 

The model describes electrons and holes via a Dirac 
equation in two dimensions using a four component 



spinor 

where the indices a and h refer to the two triangular 
sublattices of graphene and +(— ) to the two Dirac K{K') 
points. Throughout this paper, we will use the following 
representation for the Dirac matrices 

-'"^(;;)--(%o)--(o'?)^ 

stand for the Pauli matrices. 

The dynamics of the carbon background structure is 
associated with a charged scalar field (p and with a vec- 
tor field Af^ . The vacuum expectation value of the scalar 
field (ip) vanishes for pure graphene, while for doped or 
deformed graphene {(p) ^ 0. If (tp) ^ 0, the model gener- 
ates a fermion mass via spontaneous symmetry breaking. 

The lagrangian density for the gauge model reads 

+ Dt'ip^D^ip - Viip^ip) 

-{f^-'F,,, (3) 

where the fermionic covariant derivative is I?^; = 9^ + 
igA^, with a similar expression for the bosonic covariant 
derivative after replacing g by 5^, 

up to a constant Vo, and 

F^, - d^A, - d,A^ (5) 

is the usual Maxwell tensor for the vector bosonic field. 

In what concerns the electrons and its interactions with 
ip and A^, the lagrangian ([3]) is a particular case of the 
possible theories studied in^. In the following, we will 
consider only the solitonic solutions of the scalar field 
equation of motion. The lagrangian ([3]) is an example of 
fermionic dynamics. 

In ([?]), the coupling constants A4, Ag parametrize 
the self-couplings of the carbon structure and are associ- 
ated with the carbon mechanical properties like phonon 
masses, phonon self-interactions and graphene stiffness. 

The lagrangian density ([3]) is invariant under local 
C//(l) ® Uh{^) transformations. Recall that 'p'^ip can be 
identified with a dynamical fermionic mass that is a func- 
tion of space and time. Given that for doped graphene 
{•p) 7^ 0, this term generates a fermionic mass gap, which 
vanishes for the ground state of pure graphene. Besides 
the generating of a mass term for the fermions, the model 
also gives masses for the scalar and vector degrees of free- 
dom. The model allows for vortex solutions for the pure 
bosonic sector of the theory. 
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One of the goals of the present work is to investigate 
how the bosonic degrees of freedom in ^ can be excited 
via the coupling with an external electromagnetic field 
Cp. The coupling with is introduced in the usual way, 
i.e. via the minimal prescription, where a ie should 
be added to the covariant derivative, where e stands for 
a generic electric charge. 

A. Equations of Motion 

The equations of motion associated with the various 
fields are derived from C in the usual way. For fermions, 
they are given by 

[i-ff^D^ - (32 + ih2 75) (^V) = , (6) 

where I?^ = dfj, + igA^ + ie is the covariant derivative. 
The corresponding equation for the scalar field (p is 

D^Df^Lp = - g2 f tp ij - ih2 (fi ip ^5'^ 

+A4 {ip^ip)+Xe (^V)'}<P(7) 

with the covariant derivative being = 9^ + iQipA^ + 
ietpC^ and e^p being the effective electric charge of if. The 



n^T°° = ■ (\7 - igl ~ ieC^ ij + 

+ Trt TT + Dip"^ ■ Dip + V (((sV) 



gauge field equation of motion reads 
d^F"^ = -gi'^-'ij-tg^^HD''^)+ig^{D''^)^ ^. (8) 

In what concerns the coupling with the external elec- 
tromagnetic field Cp, the lagrangian density and the 
equations of motion show that, to lowest order in pertur- 
bation theory, the coupling with if is of order e^p, while A^^ 
is of the order g^pC^p or e^. Then, one expects the degrees 
of freedom associated with ip to have a high probability 
of being excited through the incidence of an electromag- 
netic wave. 



B. Hamiltonian Density 

The Hamiltonian density H can be read from the 
energy-momentum tensor 

U—ll),'ijj,cp''',:p,A^i 

It follows that 



+ + (10) 



where tt^ = d'^ip and it — d'^ip' arc the canonical momen- 
tum density fields associate with ip and (p\ respectively. 
The vectors E and B are the usual electric and magnetic 
fields associated with the gauge vector A^. 

III. SOLITON-LIKE SOLUTIONS FOR 
GRAPHENE AND NANOTUBES 

In this section we will discuss the solutions of classi- 
cal equations of motion. Let us for the moment disre- 
gard the contributions coming from fermions and gauge 
fields. The model assumes that the scalar and vector 
fields are associated with the carbon lattice background. 
The Dirac electrons are soft degrees of freedom and one 
expects that the fermion-scalar couplings are subleading 
for the dynamics of the scalar field. 

Assuming that the only dynamical field is (ys, the cor- 
responding classical configurations are obtained solving 
equation ([T])- In order to compute solitonic like solutions, 



we will consider the following ansatz 

ip{f,t) ^ (l){x-vt)e'''y (11) 
for a graphene sheet living on the xy-plane and 

^(r,^) = 0(^-«^)e™^ (12) 

for a nanotube whose symmetry axis is the z-axis. In 
order to describe the nanotube we use cylindrical coor- 
dinates r, 9 and z. 

The ansatze (|lip and (fT2)) reduce the partial differen- 
tial equation ([7]) to an ordinary differential equation in 
the variable C, where 

{X — vt, for graphene, 
(13) 
z — vt, for a nanotube. 

In terms of C, equation (I7|) becomes 

(1 - v^)(j)" = {fi^ + k^)(j) + A40^ + A60^ , (14) 
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where a prime means a derivative with respect to ^. 

The ansatz used to describe nanotubes, makes the 
equation of motion equal to ()14|) after the replacement 
of fc^ by / E? . In the following, we will explore the 
solutions of equation (|14p for graphene, but one should 



have in mind that the results for graphene have a directly 
translation into nanotubes. 

The energy associated with the solitonic waves pT|) is 
given by 



J-L^/2 ^ I 6 i 



(15) 



where and W are the dimensions of the graphene 
sheet. For a nanotube of radius R and length L, the 
energy associated with the non-linear wave is given 
by ([T5)) after the replacement of W by 27ri? and of fc^ by 
m? / . For graphene and nanotubes the lowest energy 
state is associate with u = 0, i.e. for a static solution, 
and with k = Q for graphene or n = for nanotubes. 
Furthermore, for the static solution, i.e. when w = 0, 
minimizing the energy is equivalent to solve the classical 
equation motion ([T4|) . 

From the point of view of the differential equation ([T4)) , 
the solutions describing travelling solitonic waves, i.e. for 
solutions with ?; 7^ 0, are the static solutions after the 
rescaling of the potential parameters /i^ + A:^, A4 and 
Ag by the inverse of 1 — u^. Note that for d < w 
c/300 the rescaling docs not change the nature of the 
solution derived for f = 0. However, for v > vj the 
rescaling formally changes the sign of coupling constants 
and can lead to a change on the functional form of the 
solution. Given that we aim to investigate non-linear 
waves associated with the carbon structure of graphene 
and nanotubes, one expects to be closer to the first case, 
V < Vf, and we will not explore the situation where the 
solitonic wave moves faster than the fermions v > vj. 
From now on, we will assume that v ~ except where 
stated explicitly. 

Before starting to discuss the solutions of the equation 
of motion p4|) , let us define the notation. In the follow- 
ing, the term will be included in the definition of the 
potential, i.e. we will write 



1. Kink Solution 



y((/>) = (/? + k^) <j? + ^Xi4,^ + ^K<f 



(16) 



Then, the energy becomes 



E = W dx 

rL^/2 

2W I dx V{<j)) , 



(17) 



where to write the last expression we have used the equa- 
tion of motion. 



An analytical solution for ([T4|) can be obtained when 



Afi = and 



V{4) = ^ {<t>' ■ 



(18) 



Note that a constant Vq = A402/2 was added to the orig- 
inal potential V{(j)) in order to have a positive defined 
potential energy. In terms of the original potential one 
has Ag = and + k'^ = —X^cj)^. The classical configu- 
rations are 



(/.(C) = ± 00 tanh <; 0o\/ Y (C - Co) ;> , (19) 



where is a constant of integration. The plus (minus) 
sign solution is known in the literature as kink (anti- 
kink) . For an infinite length sheet of graphene, the energy 
associated with (IT^ is given by 

E.oi^^^J^M'l^lW . (20) 



2. A Solution for a cj)^ Potential 

The classical configuration associated with the poten- 
tial energy 

V{<P) = ^c^^cf'-cflf , A>0 (21) 

was investigated by Lohe ini^. The constants in the origi- 
nal potential are given by Ag = 3A^/2, A4 = — 2A^0q and 
n"^ + k"^ = A^0o/2. The equations of motion have the 
solutions 



m = ±^Jl± tanh 



V2 



(C-Co) 



(22) 



The energy associated with any of these configurations is 
given by 



ol 



(23) 
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FIG. 1: The solutions (fT9)) and ([22} in arbitrary units. All 
constants where set to one and ("o = 0. 



for a graphcnc with infinite length and for any combina- 
tion of plus and minus signs. 

The profile of solutions (|T9|) and ((22|) can be seen in 
figure [H 



A. Energy Estimation 

The system of units used so far take ti = 1 and the 
Fermi velocity vp ~ 1. In order to convert to the 
usual system of units, we will use for the speed of light 
c = 2.998 X 10^ m/s, vp = c/300 and hvp = 0.6578 eV 
nm. To estimate the energy associated with the solitonic 
configuration, recall that in the gauge model the fcrmion 
gap is twice the fcrmion mass, i.e. 



I « 2/(/)o ^/X4 — 14/v'A4 nm, for A4 given in eV. For the 
other nonlinear wave I V2/X (j)Q = 43/A nm. 

In order to understand how the energy of the non-linear 
wave changes with W, one has to consider Eg{W). Ini^ 
Eg[W) was investigated experimentally and the authors 
conclude that for W > 16 nm, the observed gaps are well 
described by 



a 



W-W* 



(28) 



where a = 0.2 eV nm and W* = 16 nm. Taking into 
account Eg{W), it follows that the energy of the solitonic 
solutions behave as 



E 



W 



[W-W*) 



3/2 



and 



W 



{W-W*) 



2 ' 



(29) 



where the first result is for the kink solution (see section 
nil 1[) and the later for the solutions of section IIII 21 In 
both cases E vanish in the limit of infinite width, i.e. 
when graphene is recovered. This result opens the pos- 
sibility of tuning the energy of the nonlinear modes by 
chosen the nanoribbon width. 

From the point of view of exciting these nonlinear 
modes of the carbon honeycomb structure, the results 
summarized in equation p9p mean that as graphene be- 
comes wider, it will be to easier to excite the solitonic 
modes. In particular for bulk graphene these solitonic 
modes have zero energy and they must be a component 
of the graphene ground state wave function. 

As a side remark, we would like to point out that as the 
width becomes sufficiently small, the edges of the carbon 
structure, being zig-zag or armchair-shapped, should be 
taken in the computation of the solitonic energy. 



Eg ^2 



9l 



hi 



(24) 



see^ for details. For nanoribbons, a typical value for the 



gapi^ being 0.1 eV for W 
it follows that 



10 nm. For 



9l 



hi 



1, 



= 0.05 eV 



and 



0.23 eV 



1/2 



W 



10 nm 



(26) 



for the kink solution discussed in section UlI 11 where A4 
is given in eV and the width W in nm, and 



Esol 



6.7 mcV A 



W 



10 nm 



(27) 



for the soliton of section IIII 2\ recall that A is dimen- 
sionless. Note that according to ([20]) and (f24| . 
the energy associated with the solitonic configurations 
is driven by the mass gap. Furthermore, from (|19p and 
(|22p one can estimate the dimension I of the transition 
region associated with each of the solitonic solution. In- 
deed, from dyk = 0.05 eV it follows for the kink solution 



Hr = 



IV. COUPLING ELECTROMAGNETIC AND 
THE NONLINEAR WAVES 

The soliton modes are associated with a complex scalar 
field ip and, therefore, they couple with the electromag- 
netic field through and effective electric charge Cgoi- 
The Hamiltonian describing the interaction between tp 
and Cfi is given by 

/ dxdy ^ iesoi C [{df.tp'') ip - ip^ {df,p)] 

+ eL C^C^ (^V) } • (30) 

For a classical scalar field, to describe either photon 
absorption or emission, one has to evaluate the matrix 
element of Hj between the one photon state |e, q) with 
energy E-y = |g| and polarization e* and the electromag- 
netic vacuum. The decay rate of the classical nonlinear 
mode is controlled by the following matrix element, com- 
puted in the limit of long wavelengths. 



q\Hi\0) 



2'KC 



dxdy 4>{x)(l)' {x) (e-Cy) , 
(31) 
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up to a volume normalization associated with the photon 
state. In pip c stands for the velocity of light in the 
vacuum and Cy is the unit vector along y direction. 

Equation (j3ip shows that only electromagnetic waves 
polarized along the y direction, i.e. for a photon polar- 
ization along the shortest direction of the graphene sheet, 
are able to couple to the solitonic modes. For nanotubes, 
e-y is replaced by eg, the unit vector along the 9 direction, 
and the photon polarization needs an eg component in 
order to couple to the solitonic modes. 

For the kink solution fsection flll ip the integrand func- 
tion in pip is antisymmetric and the matrix element van- 
ish. In the long wavelength limit, the kink solutions do 
not couple to the photon. 

For the non-linear solution of the equations of motion 
discussed in section ITlI 21 it follows that 



(e, q\Hi\Q) = (e • e^) W. (32) 

It follows that the decay rate for photon production grows 
with the width squared and 0q, i.e. the fermionic mass 
gap squared. For the configuration ((22|) . it follows that, 
see equation ([24|). 

W,mim^ = 2^ceL- (^--e,)' ^— % (33) 

9 52 + "2 

A non- vanishing matrix element requires not only a pho- 
ton polarized along the shortest direction of the graphene 
sheet but also a solitonic wave which transport momenta 
k along the y direction. 

In what concerns the scaling of the transition rate F 
with W, combining equations ([^ and it follows 

that for ly > 16 nm 

V. CONCLUSIONS 

The gauge model for graphene suggested in^ and ex- 
tended iii^, combined with the phenomenological obser- 
vation of how the mass gap changes with the graphene 
width, predicts that solitonic waves can couple to the 
electromagnetic field, provided the photon has the right 
polarization. The decay rate is non- vanishing only when 
the photon polarization has a component along the y 
direction for graphene or along the 9 direction for a nan- 
otubc. 

From the point of view of the mass gap, the non-linear 
waves can be thought as giving rise to regions in graphene 
where the electrons behave as massless particles, those 
where 0^0, and regions where the electrons acquire a 
finite mass, regions where ~ 0o- If experimentally one 
can distinguish such regions by means of a spectroscopic 
analysis like, for example, scanning tunneling microscope 
as performed ini^ for the different graphene terminations 



or ini^ to investigate one dimensional defects, then one 
could observe directly the mass profiles and extract infor- 
mation on the solitonic waves, provided that the life time 
of these non-linear states is long enough. From these pro- 
files one can constraint the possible forms of the scalar 
potential V{(j}) and, in this way, help to better define the 
effective model. 

The problem discussed here has several energy scales 
which behave differently with W . This different scal- 
ing can, in principle, be explored to identify the soliton 
modes of section UlI 21 As discussed above. Eg oc 4%. 
Therefore, the energy associated with the soliton scales 
as 

Esoi (X E^/^W and E'^W (35) 

for the kink solution ((T^ and for (1^^ , respectively. This 
result suggests that the energy of the solitonic modes 
scales with mass gap Eg to a given power and the 
nanoribbon width. 

This different behavior can be explored, in a systematic 
investigation of the spectra of photon emission with W, 
to identify the nonlinear modes in graphene nanoribbons 
and improve our understanding of the effective theory 
discussed here. 

In what concerns the typical energy associated with 
the solitonic modes, the estimates performed at the end 
of section IIII Al suggests, for example, Esoi ^ 7 meV 
for a ^ 10 nm for the case of an effective poten- 
tial which goes as (p^. Then, the possible transitions 
involving the nonlinear modes open a window to ex- 
plore the far infrared and terahertz frequencies. Fur- 
thermore, the coupling of the nonlinear waves with the 
electromagnetic field are sensible to the photon polariza- 
tion. The nanoribbon is transparent to photons polar- 
ized along the largest length, but not to those polarized 
along the nanoribbon width. These characteristics are 
of interest for optical devices as, for example, broadband 
polarizers applied to photonics circuits for telecommuni- 
cations. Note, however, that similar properties can occur 
in association with electronic intraband transitions as de- 
scribed inii. A possible way to distinguish between the 
two mechanism is to investigate the behavior with the 
chemical potential. For the solitons, the results described 
here arc independent of the chemical potential. Dichro- 
ism with intraband transitions requires a non-vanishing 
chemical potential. 

We believe that our discussion can motivate exper- 
iments that could disentangle and identify the role of 
the carbon structure via the presence of nonlinear waves 
and/or intraband electronic transitions in grapheme 
nanoribbons. 
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